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Oj, We propose a variant formulation of Hamiltonian systems by use of variables including 

.^ • redundant degrees of freedom. We show that Hamiltonian systems can be described by 

\ extended dynamics whose master equation is the Nambu equation or its generalization. 

^^ ' Partition functions associated with the extended dynamics in many degrees of freedom 

\ systems are given. Our formulation can also be applied to Hamiltonian systems with 

I— I' first class constraints. 

^- 
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^ : 

G^l §1. Introduction 

In general, we have a choice of variables describing a physical system. In most cases, we 

^ \ choose a set of variables whose number is same as the total number of degrees of freedom of 

^^ \ the system so as to minimize the number of equations of motion. However, in some cases, 

^) ' it is quite useful to formulate the system by use of variables including redundant ones. A 

^N ' system with a gauge symmetry offers a typical example. To describe such a system keeping 

^^ . the gauge symmetry manifest, we should employ a formulation that includes redundant 

^r^ ■ variables. Although such a formulation is somewhat complicated, thanks to the symmetry, 

'^ . we can clearly understand important properties of the system such as conservation laws and 

^ \ form of interactions, and can also calculate physical quantities in a systematic way [1]. 

^ \ Therefore it is interesting to explore general features of formulations including redun- 

^ \ dant degrees of freedom. Here we are based on a principle (or brief) that physics should 

he independent of the choice of variables to describe it, and make an attempt to formulate 

Hamiltonian systems (systems of Hamiltonian dynamics) in terms of new sets of variables 

including redundant ones. What kind of dynamics describes the time evolution of the new 

variables? 

Our strategy and conjecture are as follows. Consider a Hamiltonian system described by 
a canonical doublet {q,p). Take A'^(> 3) variables {xi,--- ,X]\f) that are functions of the 
canonical doublet, and deal with them as fundamental variables to describe the system. If 
they contain redundant variables, constraints between some variables must be induced. To 
handle the constraints, Dirac formalism [2] provides a helpful perspective, where constraints 
with Lagrange multipliers are added to the original Hamiltonian. The induced constraints 
play a similar role of Hamiltonian. As for the dynamics of N variables, Nambu mechanics [3] 
is quite suggestive. In Nambu mechanics, fundamental variables form an A'^-plet, whose time 
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evolution is generated by A^ — 1 Haniiltonians according to the Nanibu equations. Combining 
the advantages of the two theories, we conjecture that there is a formulation whose master 
equation has a form of the Namhu equation or its generalization, where Hamiltonians consist 
of the original one and induced constraints. 

Nambu mechanics is a generalization of the Hamiltonian dynamics proposed by Nambu 
forty years ago [3]. In his formulation, the dynamics of an A^-plet is given by the Nambu 
equation, which is defined hy N — \ Hamiltonians and the Nambu bracket, a generalization 
of the Poisson bracket. The structure of Nambu mechanics is so elegant that many authors 
have investigated its application. However, the applications have been limited to particular 
systems such as constrained systems, superintegrable systems, and hydrodynamic systems, 
because Nambu systems (systems of Nambu mechanics) should have multiple Hamiltonians 
or conserved quantities. For example, it has been studied how Nambu mechanics can be 
embedded into constrained Hamiltonian systems [4-9] or how constrained systems can be 
described in terms of Nambu mechanics [10]. 

In this article, we show that the structure of Nambu mechanics is, in general, hidden in 
systems of Hamiltonian dynamics. That is, Hamiltonian systems can be described by Nambu 
mechanics or its generalization by means of change of variables from canonical doublets to 
multiplets. Our formulation can be generalized to many degrees of freedom systems, and the 
associated partition functions are given. We also apply our formulation to systems with first 
class constraints. Our approach would be regarded as a complementary one to the previous 
works [4-10]. 

The outline of this article is as follows. In the next section, we give a formulation of 
Hamiltonian systems using Nambu mechanics and its generalizations. As an application, 
Hamiltonian systems with first class constraints are also formulated as Nambu systems 
in section 3. In the last section, we give conclusions and discussions on the direction of 
future works. In appendix A, we derive the Nambu equation from the least action principle. 
In appendix B, we show that a Nambu system of an A^-plet can be described by Nambu 
mechanics with an A^ + r-plet (r > 1). 

§ 2. Nambu systems hidden in Hamiltonian systems 

%2.1. Review 

We begin with a brief review of Hamiltonian systems and Nambu systems [3]. A Hamiltonian 
system is a classical system described by a generalized coordinate q = q{t) and its canonical 
conjugate momentum p = p{t). These variables satisfy the Hamilton's canonical equations 
of motion, 

dq dH dp dH 

di^~df ' dt^~~d^ ' ^ ' 

where H = H{q,p) is the Hamiltonian of this system. For any functions A = A{q,p,t) and 
B = B{q,p,t), the Poisson bracket is defined by means of the 2-dimensional Jacobian, 

SA, BU ^ ^^^'^^ ^dAdB_dAdB 
d{q,p) dq dp dp dq 
In terms of the Poisson bracket, the Hamilton's canonical equation of motion for any function 
f = f{p, q) can be written by 

! = {/.//)„. (3) 
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On the other hand, a Nanibu system is a classical system described by a multiplet. As 
the most simple example, let us consider a Nambu system described by a triplet x = x{t), 
y = y{t), and z = z{t). These variables satisfy the Nambu equations 

dx _ d{Hi,H2) dy _ d{Hi,H2) dz _ d{Hi,H2) 
dt d{y,z) ' dt d{z,x) ' dt d{x,y) 

where Hi{x,y,z) and H2{x,y,z) are "Hamiltonians" of this system. For any functions A = 
A{x, y, z, t), B = B{x, y, z, t), and C = C(x, y, z, i), the Nambu bracket is defined by means 
of the 3-dimensional Jacobian, 

{A,B,CU^^^l^. (5) 

d{x,y,z) 

In terms of the Nambu bracket, the Nambu equation for any function / = f(x, y, z) can be 
written by 

f = {/,//!, F2}nb . (6) 

dt 
It is straightforward to extend the above formalism to a system described by an A^-plet Xi 
(i = 1, 2, • • • , A^). These variables satisfy the Nambu equations 

dxj _ sr-^ dHi OHn-i , . 

Clt , , OXi-^ OX'i^-^ 

where Ha = Ha{xi,X2, • • • , xj\f) (a = 1, • • • , A^ — 1) are "Hamiltonians" of this system and 
^Hi—iN-i is the A^-dimensional Levi-Civita symbol, the antisymmetric tensor with ei2---Ar = 1- 
For any functions Aa = Aa{xi,X2, • • • , X]\[,t) (a = 1, • • • , A^), the Nambu bracket is defined 
by means of the A^-dimensional Jacobian, 

SA A A \ - d{Ai,A2,--- ,An) 

{Ai,A2,---,^N}nb = -TT, T 

d{xi,X2,--- ,xn) 

y e- ■ ■ ^^...^ (8) 

. . . -, ClXi-. OXi^ C'X^^r 

2l,«2, ■■■,*« = ! 

In terms of the Nambu bracket, the Nambu equation for any function / = f(xi,X2, • • • , xjy) 
can be written by 

^^={f,Hi,H2,--- ,HN-ihB. (9) 

§ 2.2. Hidden Nambu structure 

Here let us describe a Hamiltonian system with a canonical doublet {q,p) by means of A^ 

variables rcj = Xi{q,p) (i = 1, • • • , A^). 

§ 2.2.1. Formulation. First we study the case with N = 2, for completeness. We assume 
that X = xi{q,p) and y = X2{q,p) satisfy {x,y}pB y^ 0. In this case, the equation for a 
function f{x,y) = f{q,p) is written as 

df ^ dif,H) ^ dif,H)d{x,y) ^ d{f,H) ^ 

dt d{q,p) dix,y) diq,p) d{x,y) ^"^'^^'^ ' ^ ""^ 

where H{x,y) = H{q,p). If {x,y}pB = 1, the transformation {q,p) — )■ {x,y) is the canonical 
transformation, and {x,y) are canonical variables. 
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Next we study the case with A^ = 3. We assume that variables x = xi{q,p), y = X2{q,p), 
and z = X3{q,p) satisfy at least two of the conditions {x,y}pB ^ 0, {?/, z}pb 7^ 0, and 
{z,x}pg 7^ 0. In this case, the equation for a function f{x,y,z) = f{q,p) is written as 

<i = ^iMl = ^iMl{^^ y}^^ + ^iMl{y, zU + ^^^^{Z, X}PB , (11) 

dt d{q,p) d{x,y) ' ^"^ d{y,z) ' ™ d{z,x) 

where H{x,y,z) = H{q,p). Note that q, p, and H are, in general, not uniquely determined 
as functions of rr, y, and z. 
Introducing a function G = G{x, y, z) that satisfies the conditions 



dG d{y,z) dG d{z,x) dG d{x,y) 
dx d{q,p) ' dy d{q,p) ' dz d{q,p) 

Eq. (11) is rewritten as the Nambu equation in the form of Eq. (6), 



(12) 



^ = {/,^,G},B, (13) 



where we use the formula 



d{A,B,C) _ d{A,B) dC d{A,B) dC d{A,B) dC 
d{x,y,z) d{x,y) dz d{y,z) dx d{z,x) dy 

The conditions (12) are compactly expressed as 

* j,k=l fc=l '' 

In appendix A, the Nambu equations in the form of Eq. (4) are also derived from a 
Hamiltonian system with a canonical doublet {q,p) using the least action principle. 

By use of Eq. (15), it is shown that the Poisson bracket between G{q,p) = G{x,y,z) and 
an arbitrary function u{q,p) = u{x,y,z) vanishes such that 

d{G,u,G) 



d{x,y,z) 



. (16) 



This means that G is a constant. We can eliminate the constant by redefining G, and the 
resulting G{x,y,z) = can be regarded as a constraint, which is induced by enlarging the 
phase space from {q,p) to {x,y,z). 

Here we give two comments on the induced constraint G{x,y,z) = 0. First, in case that 
dG/dz / 0, we can solve G{x,y,z) = for 2; and obtain z = z{x,y). Because the condi- 
tion dG/dz = {x,y}pp ^ also enable us to express q and p as functions of x and y, 
the expression z = z{x, y) can also be obtained by inserting q = q{x, y) and p = p{x, y) 
into z = z{q,p). Therefore the implicit form of the constraint G{x,y,z) = has an equiva- 
lent explicit form z = z{x,y), which clearly shows that z is a redundant variable in this 
case. Second, H{x, y, z) is not uniquely determined as a function of x, y, and z, that 
is, we can add a term X{x,y,z)G{x,y,z) to H{x,y,z), where X{x,y,z) is some function. 
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If a Hamiltonian H{x,y,z) satisfies H{x,y,z) = H(q,p) and Eq. (13), another Hamilto- 
nian H{x, y, z) + A(a;, y, z)G{x, y, z) also satisfies them. This is because the additional term 
X{x,y,z)G{x,y,z) always vanishes on the Nambu bracket. 

It is straightforward to extend the above formulation to the case with general A^(> 3). We 
assume that at least A^ — 1 of {xj, xj}pb (i, j = 1, • • • , N) do not vanish. In this case, the 
equation for any function /(xi, • • • ,X]y) = f{q,p) is written as 

df_dif,H)_l^dif^ 



dt 



d{q,p) 2 f-^ d{xi,Xj) 



where H{xi, • • • , xn) = H{q,p). 
Introducing functions Gb = Gb{xi, ■ ■ ■ , xn) {b = 1,- ■ ■ ,N — 2) that satisfy the conditions 

1 v^ d{Gi,--- ,Gn^2) _ f , , s 

(iV-2)! . ^ ^i^i^i^-i- 5(xi3,--- ,Xi„) -i^^^-^^^tpB , U»j 

Eq. (17) is rewritten as the Nambu equation in the form of Eq. (9), 

-TT = {fyH,Gi,- ■ ■ ,Gn-2}nb , (19) 

where we use the formula concerning Jacobians, 

d{AuA2,--- ,An) _ 1 ^ d{Ai,A2)d{A3,--- ,ijv) ( . 

dix,,X2,--- ,xn) 2(iV-2)!..^. /^^'^^^■■■*"9(x,„x,J9(a;,3,---,x,J ■ ^ ' 

By use of Eq. (18), it is shown that the Poisson bracket between any of A^ — 2 functions 
Gb{Q,p) = Gb{xi,X2, ■ ■ ■ , xn) and an arbitrary function u{q,p) = u{xi,X2, • • • , xn) vanishes 
such that 



fr< \ - ^ sr^ d{Gb,u) ^ . 

\<-^6) 'U./PB — 7J / ^ TTT \lXi11Xi2JpB 

1 ^ d{Gb,u) djGi,--- ,Gn-2) 

2{N-2)l ^ (^^ '^^^^'^■■■^"9(x,„x,J 9(x,3,...,x.J 

^ d{Gb,u,G\,--- ,Gn-2) ^ PI /2i ^ 

5(xi,X2,X3,--- ,xn) 

Hence Gb are constants. We can eliminate the constants by redefining G^, and the resulting 
Gfe(xi,X2,--- ,X7v) = can be regarded as induced constraints, which are associated with 
enlarging the phase space from {q,p) to (3;i,3;2, • • • ,xn)- 

In this way, Hamiltonian systems can be formulated as Nambu systems by use of N variables 
Xi = Xi{q,p) (i = 1, 2, • • • , N). The variables form an N-plet, and the N — 1 Hamiltonians 
are given by the original Hamiltonian H{xi,X2, • • • ,xn) = H{q,p) and induced constraints 
G'fe(xi,X2, • • • ,X7v) = (6 = 1,--- ,iV — 2). Note that H{xi,X2,- • • ,xn) is not uniquely 
determined, because of the freedom to add a term ^^ Afe(xi, X2, • • • , XAr)G'fe(xi, X2, • • • , xn) 
to i/(xi, X2, • • • , xtv). Here Afe(xi, X2, • • • , x^v) are some functions. 
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§ 2.2.2. Examples. Here we present two simple examples to show how induced constraints 

are obtained for given multiplets. 

(a) A^ = 3 

Consider composite variables, 



^ '" /) , y = -Aq^+p') , z=\qp, (22) 



which satisfy the following relations, 

{x,y}m = z, {y,z}pB = x, {z,x}pB = -2/. (23) 

Then the conditions (Eq. (12)) become 

dG__ dG__ dG_ 

dx dy ^ dz 

and G is obtained by 

G=^{x^-y^ + z^)+C , (25) 

where C is a constant. Redefining G as G — C, we obtain the induced constraint G{x, y, z) = 
G{q,p)=0. 

(b) iV = 4 

Consider variables including composite ones, 

xi= q , X2=p , X3 = X3{q,p) , X4 = Xi{q,p) , (26) 

which satisfy the following relations, 

(fT'o (J'T A 

{xi,X2}pB = l, {a;i,X3}pB = -^ , {a;i,X4}pB = -^ , 
(9x3 dxi 

|a:2,X3JpB = --^ , {X2,X4|pB = — -K- ■, 

r . _ dx3_dx4_ _ dx3_dxi 

dq dp dp dq 



Then the conditions (Eq. (18)) become 



E dGi dG2 _ . , , . 

^iiiii-iiA p. Pi — X^iiiXiiiPB 1 l^oJ 



and Gi and G2 are given by 

Gi = X3 - 2;3(j;i,X2) + Gi , G2 = X4 - X4(xi,X2) + G2 . (29) 

where Gi and G2 are constants. By means of redefinitions of Gi and G2 to elimi- 
nate the constants, we obtain the induced constraints Gi(xi,X2,X3,X4) = Gi{q,p) = and 
G2(xi,X2,X3,X4) = G2(g,p) = 0. 

§ 2. 3. Many degrees of freedom systems 

Let us extend our formulation to Hamiltonian systems with many degrees of freedom. 
Consider a Hamiltonian system described by n sets of canonical doublets iq(k),P(k)) (k= 
1, 2, • • • , n). As is the case with n = 1 given in § 2.2, hidden Nambu structure can be also 
found in this system. Here we present A^ = 3 case, that is, the case with n sets of triplets 
Xi(]^) = 2;i(k)((7{k),P{k)) {i = 1,2,3). Generalization to A^(> 3) cases is straightforward. 
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i^ 2.3.1. Dynamics. In this system, the Poisson bracket of A and B is defined by 

_^ / dA dB dA dB \ 

and the Hamilton's equation of motion for any function / = f{q(i),P(i), • • • , Q{n),Pin)) can be 
written as 

f = {/,^}pB, (31) 

where H = i7(q'(i),P(i), • • • ,(?(„), P(„)) is the Hamiltonian of the system. On the other hand, 
the Nambu bracket oi A, B, and C is defined by 

{A,B,CU^±-^^^1^^, (32) 

where a:(k) = xi(k), y(k) = a;2(k), and 2;(k) = X3(k). Then the Nambu equation for any function 
/ = /(a;(i),y(i),--- ,2;(„)) can be written as 

^ = {f,H,GU. (33) 

Here ^ = ^(x(i),y(i,, • • • , Z(„)) = iJ(g(i),P(i), ■ ■ ■ , g'(„),P(„)) is the Hamiltonian and G = 
G(a;(i), 2/(1), • • • , 2;(„)) = ^j, G(k)(x(]j), t/j^), Z(k)) is a sum of the induced constraints that satisfy 
the conditions 

^G'(k) _ a(j/(k), ^(k)) 5G'(k) _ ^(^(k);3;(k)) oG^t^^ _ t/(x(it),|/(k)) , ,^ 

9X(k) 5(g(k),P(k)) ' %k) 9(g(k),P(k)) ' 5Z(k) 5(g(k),P(k)) 

Note that the induced constraints are defined so as to be zero, G'(i,)(x(k),y(k);-2{k)) = 
G'(k)(^(k)5P(k)) = 0, and the Hamiltonian is not uniquely determined because of the freedom 
to add a linear combination of G(k) to H. 
The 3n variables x^t^^ satisfy the relations 

|2;ji{ki)5 2;J2{k2)) 2;i3(k3)}NB = £iiJ2«3 for ki = k2 = ks , (35) 

{2;ji{ki),a;i2{k2),2;i3(k3)}NB = for otherwise. (36) 

The first type of relations (Eq. (35)) are invariant under the time evolution (Eq. (33)) 
irrespective of the form of H. To be more specific, for infinitesimal transformations ^^(k) -^ 

l^(k)'2/(k)l ^(k)}NB = 1 (37) 

hold. We can also show an important relation, 

(^{^(1) 1 ^(1) 1 ^(1) i • • • ; ^(n) 1 j/(n) ; ^(n) ) ^ ^ ,gg^ 

C'(3^(l) 5 ?/(l) 5 ^(1) ) ■ ■ ■ ; ^(n) 5 2/{n) ) ^(n) ) 

which guarantees the Liouville theorem, the conservation law of the phase space volume 
under the time development. On the other hand, the second type of relations (Eq. (36)) do 
not always hold, unless there is no interaction between the n subsystems, i.e., H has the 
form such as #(x(i), y(i), • • • ,Z(„-,) = Ek-^(k)(2;{k),2/(k), %))• 
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§ 2. 3. 2. Partition functions. It is well-known that the partition function Zh for a canonical 
ensemble of the Hamiltonian system ((7(i),P(i), • • • ,(l{n)-,P(n)) is defined by 



Zh = 



J|dg(k)dp(k)e ^" 



(39) 



k=i 



where (3 = \/{kBT) is the inverse temperature made of the Boltzmann constant ks and the 
temperature T. Here we study the partition function Zn for an ensemble of the Nambu 
system (a;(i), y(i), • • • , 2:(„)) hidden in the Hamiltonian system. 

First let us conjecture the form of Z^ on physical grounds. Since H = H, Z^ must contain 
the "Boltzmann weight" such as e~" . The other Hamiltonian G is a sum of constraints 
G(k) (a;(k) , 2/(k) ) -^ck) ) = G'(k) (^(k) , P(k) ) = 0, and therefore there should be the delta functions such 
as 5(G(k)) in Zn. Furthermore, Zn must contain the volume element nk=i dx(]^)dy(]^)dz(]^) from 
the Liouville theorem. 

On the basis of above observations, it is expected that Zn should have the form such that 



Zn 



JJda:(k)d2/(k)d2:(k)(5(G'(k))e '^^ 



k=i 

n 



JJ(iX(k)(iy(k)d2;, 



(k) 



k=l 



d7, 



27r 



Hg-Z^^-^^kTc^Gck) 



(40) 



(41) 



We can derive Zh (Eq. (39)) from this expression of Z^. For example, let us consider the 
case that 9G'(k)/5z(k) / 0. We assume that there are A'k solutions of G'(k) = 0, z^^ (ok = 
1, 2, • • • , A'k), and all of them satisfy the conditions (Eq. (34)). Then using the formula of 
the delta function and the change of variables, Eq. (40) becomes 



Zn 



A^k 



_[_[ t?X(k)(iy(k)ii2;(k) 2_^ "^1%) ~ ^(k) \X{yi),y('n))) 

ak=l 

A^k 

H (ix(k)(iy(k)(i2;(k) 2_^ "^l^) ~ ^(k) {^WiVo^))) 

k=l ak=l 



k=l 
n 



dG 



(k) 



dz. 



(k) 



-I3H 



(j\^{k)-,y(k)) 



d{qw,Pw) 



-PH 



M 



jj-j n rf'?(k)dP(k)e-^^ = AAZh 

k=l 



(42) 



where J\f = nk=i -^k is a constant normalization factor. This factor is irrelevant to the 
evaluation of physical quantities. 

It is natural to require that Zn should agree with Zh (up to some normalization factor), 
because we just describe the same physical system using different formulations. It should be 
noted here that either expression of Zn (Eq. (40) or Eq. (41)) is different from that proposed 
in Ref. [3]. It would come from the fact that the Nambu mechanics considered here is an 
effective one induced due to the redundancy of the variables. 

Finally we just give the result for the case of general A^(> 3). The possible form of the 
partition function is given by 



Zn 



Y\ dxn^)dx2(k) ■ ■ ■ da;Ar(k)5(G'i(k))(^(G'2(k)) • • • 6{GN-2(k))e~ 



I3H 



(43) 



k=l 



where Gb{k) =0 (6= 1,2,--- ,N — 2) are induced constraints. This expression should agree 
with Zh (Eq. (39)) up to some constant normalization factor. 
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^2.4- Generalized Nambu equations 

We generalize our formulation to include a specific case that all multiplets share some 

variables. In such a case, a generalization of the Nambu equation would be required. 

Let us describe a Hamiltonian system with n sets of canonical doublets (9(k))P(k)) (k = 
1, • • • , n) using 2n + ni variables we (£ = !,••• , 2n + ni). We classify the variables wi into 
two groups, Xa {a = !,■■■ , 2n) and Zs (s = !,••• ,m) where Xa are assumed to satisfy 
det{xa,Xb}pB 7^ 0. Note that the classification of variables is not unique. 

First we study the case with tti = for completeness. In this case, the equation for any 
function f{xi, ■ ■ ■ , j;2n.) = fiqw,P(i), ■ ■ ■ ,9(n),P(„)) can be written as 

df _y d{f,H) _ly y d{f,H) d{xa,xt) y -^^EMl (44) 

where H = H{xi,--- , X2„,) = H{q<^^ , p(i) , • • • , q^„^ , p(„) ) and gab is defined by 

gab{xi,--- ,X2n) = 5afe(g(l),P(l) , ' ' ' ,5(n),P(n)) = o T] ^7 ^^^ 7 = l^{Xa,Xb}pB ■ (45) 

^ k=l '^W(k),P(k)J ^ 

The ^ab plays a role of the metric tensor, because it transforms under a change of variables 

Xa — > x'a as follows, 

g'abix'i,--- ,4n)= 2-. 'dfdx'^'"^''^^''" '^^'''' ■ ^'^^'^ 

c,d=l 

In case that ^aj, depends on Xa, neither the transformation (9(i),P{i), • • • ,q(n),P(n))^- 
(xi, • • • , X2n) nor the time evolution of Xa is the canonical transformation. The latter means 
that the Liouville theorem in general does not hold for the dynamics of Xa ■ This fact reminds 
us superiority of canonical variables. 

Now let us proceed to the case with m > 1. The equation for a function 
f(wi, ■ ■ ■ ,W2n+m) = f{q(i),P(i), ' ' ' ,q(n),P{n)) Can be written as 

2n 



di _ 1 y^ d{f,H) 

a, 0=1 



2n m 



where H = H{wi,--- ,ty2n+m) = -f^(Q'(i),P(i), • ' ' ,Q(n),P{n))- Introducing functions G^ {s = 
1, • • • , m) and g^ that satisfy the following relations, 

-{Xa,XbU = 9ab diz„...,Z^) ' ^"'^ 

2n ~ ^ ~ ~ ~ 

Ir 1 Sr^ ~{m)d[Gl,- ■ ■ ,Gs~l,Gs,Gs+l,- ■ ■ ,Gm) /.„N 

2 f^ d{zi,--- ,Zs-l,Xb,Zs+l,--- ,Zm) 

2^ ~ ~~~ ~~~ ~ 

r T V^ ~(m) <9(Gi, • • • ,Gs~l,Gs,Gs+l,- ' ' ,Gt~l,Gt,Gt+l,--- ,Gm) fr-n\ 

{Zs,Zt}pB = > g'ab W-, ^ , 50) 

ab^l 0[Zi,--- ,Zs-l,Xa,Zs+l,--- ,Zt~l,Xb,Zt+l,--- ,Zm) 
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where s < t, Eq. (47) can be rewritten by 



df _ sr^ ~(m) d{f, H,Gi,- ■ ■ , Gm) ,^^x 



using a formula concerning Jacobians. 

By use of Eqs. (48)-(50), it is shown that the Poisson bracket between any of 
TTi functions Gs{q(i),P(i),- ' ' ,Q{,i),P(n)) = Gsiwi,- " ,W2n+m) and an arbitrary function 
uiqii),P(i), ■ ■ ■ ,q(n),P{n)) = u{wi, ■ ■ ■ ,W2n+m) vanishes such that 



2n 






2n m 



X; 5ir^ diGs,u,G,,.--,G^) ^ Q _ (52) 

a, 0=1 



Hence Gs{q(i),P(i), • • • ,q(n),P(n)) are constants and if necessary we can define Gs = Gs = 
by shifting constants. We refer to Eq. (51) as the generahzed Nambu equation. Note that 
the Liouville theorem does not hold in general for the dynamics given by this equation. 
This unfavorable property comes from two reasons: Eq. (51) has Xa-dependent g^ and 
multiplets in Eq. (51) share common variables Zg. The latter means that it is difficult to 
define appropriate phase space volume. 

One of the non- vanishing components of g^ can be set to g by redefinition of constraints 
Gg- For example, in case that n = 1, we can set g]^ = ^ (and 521^ ~ "2) ^^ redefining Gg, 
and Eq. (51) reduces to the Nambu equation (Eq. (19)) with N = 2 + m. 

Finally we consider the case that variables Xa and Zg are further classified into 
M "irreducible" sets, {x^V , zf,^ } © {x^^^ , zfj }®---®{x[f,, z^^t^ } , where a' = !,■■■ M 
(Si=i '^^ ~ '^) ^'^'i s* = 1, • • • , m* (X]i=i "^ — '^)- Here "irreducible" means that the Poisson 
bracket between any two elements that belong to different sets vanishes, i.e., {x^'- , x^^] }pb = 0, 
{x^'i ,z^''j'}pB = 0, and {z^l', z^j'jpe = for i^ j. Note that this classification is not unique, 
either. The equation of motion for any function /(t^i, • • • ■,W2n+m) can be expressed in the 
form of the generalized Nambu equation, 



^^Sr Sr --(m-) d{f,H,G^ r-- ,G^^) 

"^ i=la',fe'=l '^Va^a"%'^l ' '^m'i 
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Here G^- and g^JlJ should satisfy the following conditions, 
lr„W ^{')T _ ~{m') diG; ,■■■ ,G;^0 

^ ^1^1 ' ■ ■ ■ ' vo 

b- = l C'l^l )■■■ ,^s'-l'^fe''^sHl''" '^m'^ 



a\b^=l 






where s* < t*. We refer to the systems where the master equations are given by Eq. (51) or 
Eq. (53) as generalized Nambu systems. 

§ 3. Nambu systems in constrained Hamiltonian systems 

%3.1. Subject 

In the previous section, we found that a Hamiltonian system can be formulated as a Nambu 
system with multiplets including composite variables of q and p. The main feature of our 
formulation is the existence of induced constraints that are required just for the consis- 
tency between the variables. Together with the Hamiltonian of the original system, the 
induced constraints serve as Hamiltonians of the Nambu system. Therefore it is intriguing 
to study how constrained Hamiltonian systems, systems with physical constraints, are cast 
into Nambu systems in our formulation. 

Relations between Nambu systems and constrained Hamiltonian systems have been inves- 
tigated by many authors [4-10]. To clarify the difference between previous works and our 
approach, here we give a brief summary of the results obtained so far. In most works, Nambu 
systems are treated as the original systems, and studies have been carried out to find appro- 
priate constrained Hamiltonian systems which the Nambu systems can be embedded into 
[4-9]. Specifically, it has been shown that Nambu equations (Eq. (4)) are compatible with 
the following equations, 

m = H^^^, (57) 

V^^i^^ = 0. (58) 

Here pi {i = 1, 2, 3) are the canonical conjugate momenta defined by pi = dL/dxi with the 
Lagrangian 

1=1 
Equation (58) can be derived as the Euler-Lagrange equation from this Lagrangian, and Eq. 
(57) leads to the relations (pi ^ Pi — HidH2/dxi = 0, which can be regarded as constraints. In 
this way, Nambu systems can be interpreted as Hamiltonian systems with specific constraints. 
On the other hand, it has been studied whether constrained systems can be described as 
Nambu systems or not. Specifically, it has been shown that constrained Hamiltonian systems 
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can be formulated in terms of (a generalized form of) Nambu mechanics by introducing an 
extra phase-space variable [10]. For a system with canonical variables {qk,Pk) (^ = 1) ' ' ' ^n) 
and m first class constraints 0/((?i, • ' ' ^Pn) = 0, the equations of motions are given by 

dqk dH v-^ / dXi , , dchi \ , , 

dt dpk ^ V 9pk opkj 

dpk dH ^ f dXi d(pi \ 

dt dqk ^ V 9qk dqk J 

where A; are Lagrange multipliers. Equations (60) and (61) are derived from (a generalized 
form of) Nambu equation 

df _^ d{f,H^,H2) 

dt ^ d{qk,Pk,r) ' ^ ^ 

where / = f{qi, ■ ■ ■ ,Pn), ^ is an extra phase-space variable, and Hamiltonians are defined 



as 



m 



Hi = H-r , H2 = r+y Xi(pi . (63) 



1=1 



The equation for r is given by 



§ = - E (^'i'^'' ^}™ + U\u HU) = -Yl ^'^' (64) 

1=1 1=1 

where the last equality holds after imposing constraints. Requiring the extra variable r 

should decouple from the dynamics, i.e., dr/dt = 0, we obtain d(j)i/dt = 0. 

Our approach differs from these previous works. Our starting point is not Nambu systems 

but Hamiltonian systems with constraints, and we do not introduce extra variables but use 

redundant variables. 

§ 3.2. Nambu structure in constrained Hamiltonian systems 

Here we demonstrate that systems with first class constraints can be formulated as Nambu 

systems or generalized Nambu systems, without introducing extra degrees of freedom. 

As a warm-up, we consider a system of two canonical doublets {qi,pi) and {q2,P2) with 
one first class constraint (f){qi,pi, q2,P2) = that is time independent: d(j)/dt = {(j), /i}pb = 0. 
Here h = h{qi,pi,q2,P2) is the Hamiltonian of this system. The constraint (j) is associated 
with gauge degrees of freedom, and an auxiliary condition x{q1jPi^Q2,P2) = such that 
{(/>, x}pB / should be imposed to fix the freedom. 

By an appropriate canonical transformation {qi,pi, q2,P2) -^ iQi,Pi,Q2,P2), we can elim- 
inate one of the canonical variables. Here we show the case that P2 is eliminated as 
follows,^ 

P2 = x{qi,pi,q2,P2) = . (65) 



^ The canonical transformation generated by the generator G ^ Xcj) is the gauge transformation, 
and the infinitesimal one is given by 6qr = {dG/dpr)S^ and 5pr = —{dG/dqr)5^ [r = 1,2). Here A 
is an arbitrary function of canonical variables and 5S, is an infinitesimal parameter. If we take A — 
{P2 — x)/((^0/^'Z2)C) using a finite parameter ^, p2 is transformed into P2 ~ X under the constraint 
= 0. In the same way, one of the canonical variables can be eliminated by an appropriate canonical 
transformation. It depends on physical systems which variable should be eliminated. 
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The new Hamiltonian K is given by K{Qi,Pi,Q2, P2) = h{qi,pi,q2,P2), and the orig- 
inal constraint cj) is transformed as ^{Qi,Pi,Q2, P2) = (I){q1tP1tQ2,P2)- From {(J),x}pb = 
d^/dQ2 7^ 0, the constraint <l> = can be solved by Q2 to give Q2 = Q2{Qi,Pi)- Then we 
obtain a constraint ^ = Q2 — Q2{Qi, Pi) =0, which is equivalent to the original constraint 
(p = 0. 

If we consider a system described by the variables (Qi, Pi, Q2) with the constraint ^ = 0, 
it is easy to show that the equation of motion for any function / = f{Qi, Pi,Q2) can be 
written in the form of the Nambu equation, 

dl ^ d{f,H,^) 

dt d{Qi,Pi,Q2) ' ^ ' 

where H = H{Qi,Pi,Q2) = K{Qi,Pi,Q2,P2 = 0) is the Hamiltonian. In fact, for f = Qi 
and / = Pi, usual Hamilton's canonical equations of motion 

dQi ^ dH{Qi,P^,Q2{Qi,Pi)) m ^ dH{Qi,Pi,Q2{Qi,Pi)) 
dt dPi ' dt dQi ' ^ ^ 

are derived from Eq. (66), and for / = ^, we obtain time independence of the constraint, 
d^ /dt = 0. On the other hand, for f = Q2, the following equation is derived, 

dQ2 ^ 9(g(Qi,Pi,Q2),^) ^ d{Q2iQi,Pi),H{Qi,P,,Q2iQi,Pi))) 
dt 9(gi,Pi) 9(Qi,Pi) ■ ^ ' 

Using d^ /dt = and Eq. (68), we obtain the Hamilton's equation of motion for Q2{Qi^ Pi); 

dQ2{Qi,Pi) ^ d{Q2{Qi,Pi),H{Qi,Pi,Q2{Qi,Pi))) 

dt 9(Qi,Pi) • ^ > 

By referring the results in § 2.2, let us formulate this system by means of composite triplet 
X = X{Qi,Pi), Y = y(Qi,Pi), and Z = Z{Qi,Pi) imposing a constraint G{X,Y,Z) that 
is equivalent to the original constraint (j). We assume that d{X,Y)/d{Qi,Pi) 7^ 0, i.e., Z is 
a redundant variable such that Z = Z{X,Y). If the variables satisfy the conditions 

dG _ diX, Y) dG _ d{Y, Z) dG _ d{Z, X) 



dz d{Qi,Pi) dx d{Qi,Pi) dY d{Qi,Pi) 

then the time evolution of any function / = f{X, Y, Z) is given by the Nambu equation, 

df dif,H,G) 



(70) 



dt d{X,Y,Z) ' 



(71) 



where H is equal to the original Hamiltonian, H{X,Y,Z) = H(Qi,Pi,Q2). We can define 
various types of Nambu systems depending on the choice of variables and the constraint. 

Here we present two simple examples. First if we choose X = Qi,Y = Pi, Z = Q2{Qi, Pi), 
and G{X,Y,Z) = ^(Qi,Pi,Q2), Eq. (71) clearly holds from Eq. (66). Next let us choose 
Y = Pi, Z = Q2iQi,Pi), and G{X, Y, Z) = <^iQi,Pi,Q2)- In this case, if the variable X is 
given by 

then the variables satisfy the conditions (Eq. (70)), and the system is described as a Nambu 
system. 
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It is straightforward to extend the above "warm-up" discussion to many degrees of freedom 
systems. Consider a system of n sets of canonical doublets {qk,Pk) {k = I,- ' ' )^) with m 
kinds of first class constraints (f>s{qi,pi,- ■ ■ , qn,Pn) = (s = 1, • • • , m). The Hamiltonian of 
this system is given by h{qi,pi, ■ ■ ■ , qn,Pn)- To fix the gauge degrees of freedom, m kinds of 
auxiliary conditions Xtiqi,Pi, ■ ■ ■ ,qn,Pn) = (t = 1, • • • ,m) that satisfy det{(ps,Xt}pB / 
should be imposed. 

By an appropriate canonical transformation {qi,pi,- " ^QniPn) ^- {Qi,Pir " > Qn,Pn), we 
can eliminate some of the canonical variables. Here we show the case that Pn-m+t are 
eliminated as follows, 

Pn-m+t = Xt{qi,Pl,--- ,qn,Pn) = . (73) 

The new Hamiltonian is given hy K = K(Qi, Pi, • • • , Qn, Pn) = h{qi,pi, ■ ■ ■ , qn,Pn), and the 
original constraints (ps are transformed as $s(Qi,Pi,--- ,Qn,Pn) = 4's{qiiPi,'" iQn,Pn)- 
Prom det{(j)s,Xt}pB = det{d^s/dQn-m+t) ¥" ^1 ^^^ constraints <&s = can be solved by 

^n—m+t to give Lln—m+t — njn— m+tvn;!) ^li ' ' ' i ^n—mi ^n—m)- 

By referring the results in § 2.4, let us formulate this system by means of 
composite variables Xa = Xa{Qi,Pi, • • • , Qn-m, Pn-m) (a = 1, • • • ,2n — 2m) and Zg = 
Zs{Qi,Pi,--- ,Qn-m,Pn-m), imposing the constraints GsiXi,- ■ ■ ,X2n-2m,Zi,- ■ ■ ,2^) 
that are equivalent to the original constraints cpg- We assume that det{Xa, ^fe}pB ¥" 0) where 
the Poisson bracket is defined by 

This means that Zg are redundant variables such that Zg = Zs{Xi, • • • , X2n-2m)- 
If Gg and g^ satisfy the following relations, 

Ij-v^ 7 1' — _ \^ ~{m)d{Gi,- ■ ■ ^Gs-l,Gs,Gs+l,- ' ' -.Gm) ,„„^ 

-i^a,^s J-PB — 2-^ 9ab "afy 7 V 7 7~\ ' ^ ' 

h^l O'l^Zi, • • • ,Zs_l, Afe,Zs+i, • • • ,Zmj 

2n—2m i^//^ i^ r^ r^ r^ r^ r^ r^ \ 

i^S'-^tJPB — /, 9ab at 7 7 Y 7 7 v 7 7 \ A' ' ) 

T^i CvZl,--- , Zs-l, Aa, Z^+l, • • • , Zt_l, Afe, Zi+i, • • • , ZmJ 

a, 0=1 



5 ^m) 



where s < t, then the time evolution of any function / = f{Xi, • • • , X2n-2m, ^ir " ) ^i 
can be written by 



^ _ v^ -(m) d{f,H,Gi, ■ ■ ■ ,Gm) , . 

dt ~ ^^ ^"' 5(X„,X6,Zi,---,Z^) ' ^ ^ 

where if is the Hamiltonian, 

-"- V'^l ' ^1 ) ' ' ' 5 ^n—rm ^n~mj ^jn— m+1 ) -'n— m+1 — U, • • • , t^n, in — "j ■ \ ' >^j 

We can define various types of Nambu systems depending on the choice of variables and 
the constraint. Por example, in case that we take Xn-m+a = Pa (a = 1, • • • ,n — m), Zg = 
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Qn-m+sj and Gs = ^s, if the variables Xa are functions of Qa that are given by 

^a — ^ay^a) — ■^ Ha n-m+a a( y ... y \ "^Jo 



2 / 9rLm^.J,77^^^^^^''Tr^dQ^ . (80) 



then the variables satisfy the conditions (Eqs. (75)-(77)), and the system is described as a 
Nambu system. 

In this way, Hamiltonian systems with first class constraints can he described as Nambu 
systems where master equations are Nambu equations or generalized ones. It is straight- 
forward to formulate constrained Hamiltonian systems as Nambu systems where both of 
original constraints and induced ones serve as Hamiltonians. Such systems can be realized 
by introducing many more redundant variables. 

§ 3.3. Example: relativistic free particle 

As an example, we consider a relativistic particle moving freely in 4-dimensional Minkowski 
space. The motion is expressed by the space-time 4-vector q^^ = q^{T) and corresponding 
canonical momenta p^ = Pfj,{T) {fj, = 0, 1, 2, 3), where r is the proper time. Here we use the 
metric tensor r]^^ = diag(l, —1, —1, —1). The system has four sets of canonical pairs {q^,Pfj) 
and one first class constraint 

cj) = pf'p^ - m^c^ = , (81) 

where m is the mass of the particle, c is the speed of light, and the Einstein's summation con- 
vention is used. We impose an auxiliary condition x = 5*^ — cr = to fix the gauge freedom. 
Performing a canonical transformation 

Po^ Po= Po, Pi^ Pi=Pi , (82) 

where i = 1,2, 3, we can eliminate Q^, and the system is described by three sets of canonical 
pairs (Q*, Pi) with the new Hamiltonian K = —cPq. The original constraint (p is transformed 
as 

(f)^<^ = P^'P^ - rr?(? = , (83) 

which has an equivalent expression. 



^ = Po + yP^TmV = , (84) 

where P = ^^ P]^ . Then the Hamilton's equations of motion for Q* and Pi are given by 

"«■ ^'^ '^^" (85) 



2„2 



ra'-c 



dt dPi /p2^ 

Using the results in § 3.2, let us construct Nambu systems that are equivalent to this 
system. The target equation is Eq. (78) with a,h = 1,- ■ ■ ,6 and m = 1. Here we present 
three types of Nambu systems. In each case the Hamiltonian is given hy H = K. 
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(a) First we consider the simplest construction, 

Xi = Q , -^i+3 = Yi = Pi , Z = Pq , 



G = q, = Z + Vy^TmV , (87) 

where Y = '^^Y^. Prom Eq. (75) we obtain 

and Eq. (78) becomes 

df ' dif,H,G) 



dt 



E 



-d(Xi,Y,,Z) 



This equation reduces to the Hamilton's equations of motion (Eqs. (85)-(86)) and the energy 
conservation law . 



Xi — Q , Xj_|_3 — Yi — Pi , Z — Pq , 



(b) Next we consider a slightly different construction, 

3 = Yi = f 

G = ^ = Z-Y'^- rr?(? . (90) 

In this case, g\2 has the same form as Eq. (88), because in both cases of (a) and (b), 
{Xj,Xj_|_3}pg = 1 and dG/dZ = 1 holds. Therefore the resulting equation is same as Eq. 
(89). 

(c) Finally we consider a case that dG/dZ ^ 1, 

Xi = 2PqQ^ , Xj+3 = Yi = Pi , Z = Pq , 

G = ^ = Z'^ -Y'^ - rri^c^ . (91) 

From Eq. (75) each component of the factor ^^^ is determined as follows, 

~(i) _ ^jYj - XjYj 
9ij - 2^2 

~(l) ^X ~W - ^x 

9ii - 2"^+3,l ; 9ii - -2°l;i+3 ' 

t2 = , (92) 

where i,j = 1, 2, 3 and /, m = 4, 5, 6. Although this is different from Eq. (88), we obtain the 
same equation as Eq. (89) again. This is because d{f,H,G)/d{Xi,Xj,Z) = holds in this 
case. 

§ 4. Conclusions and future works 

We have given a variant formulation of Hamiltonian systems in terms of variables including 
redundant degrees of freedom. By use of a non-canonical transformation that enlarges phase 
space from {q,p) to (xi,X2, • • • ,X]\f), we can reveal the Nambu mechanical structure hidden 
in a Hamiltonian system. Hamiltonians required for the Nambu mechanical description are 
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given by the Haniiltonian of the original system and constraints induced due to the consis- 
tency between the variables. Our formulation can be extended to many degrees of freedom 
systems and systems with first class constraints. Generalized forms of Nambu equation (Eqs. 
(51) and (53)) are required in some cases. Our approach to constrained systems is different 
from the preceding works [4-10], that is, we treat Nambu mechanics as effective mechanics, 
and we introduce not extra degrees of freedom but redundant degrees of freedom. 

Our formulation would not just be a change of description, but give a new insight for sta- 
tistical or quantum mechanical treatment of Hamiltonian systems. For example, the Nambu 
equation (Eq. (19)) could give a basis of a novel quantization scheme for a Hamiltonian 
system. In the present work, the constraints (Gi, (^2, • • • , Gn) are unphysical ones and they 
all are set to zero. However, if we give them some appropriate values, the equation (19) could 
provide semi-classical equations for quantum-mechanical expectation values [11, 12]. Non- 
vanishing Gb come from quantum fluctuations, for example, if we take x = {q), y = {p), and 
z = (g^), then G = z — x^ has non-zero value in general. The same argument would holds for 
statistical-mechanical expectation values. Therefore we expect that the Nambu equation (Eq. 
(19)) with non-vanishing Gb could be a master equation for statistical or quantum mechanics 
of Hamiltonian systems. More detailed studies will be presented in future publication, and 
they might provide important clues on handling statistical or quantum mechanics of Nambu 
systems [13-23]. 
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A. Derivation of Nambu equation from the least action principle 

Let us derive the Nambu equation (Eq. (19)) from a Hamiltonian system using the least 
action principle. Here we show the case with A^ = 3, where the Nambu equations are given 
in the form of Eq. (4) . 
Our starting point is the action integral such that 

S . / (,| _ «„.„) ., . / (,(., , .,|,(., , .) _ ^(., , .)) * , (A,) 

where x = x{q,p), y = y{q,p), z = z{q,p), and we assume that the Hamiltonian H can be 
expressed by 

H{q,p)=H{x,y,z) . (A2) 

As mentioned in §2.2, in general, g, p, and H are not uniquely determined as functions of 
X, y, and z. Our goal is to obtain the equations of motion that hold independently of the 
expressions of q and p. 

Let us regard x = x{t), y = y{t), and z = z{t) as independent variables, and consider their 
variation x ^- x + 5x, y ^ y + 6y, and z — ?• z + 6z. Integrating by parts and ignoring the 

17/20 



surface terms, the variation of S can be written by 

^Sx + —6 + —6z] — - 
dx dy dz J dt 

( dH ^ dH ^ dH \ 
\ ox oy oz I 

(dp dq dq dp dH \ 
dx dt dx dt dx j 

I dp dq dq dp dH \ 
\ dz dt dz dt dz J 



dq^^ dq^ dq^\ dp 
dx dy dz J dt 



dt 



dp dq dq dp dH \ 
dy dt dy dt dy I 



dt . 



(A3) 



Imposing the least action principle 6S = 0, we obtain the equations of motion 



d{q,p) dy d{q,p) dz 

d{x,y) dt d{z,x) dt 

d{q,p) dz d{q,p) dx 

d{y,z) dt d{x,y) dt 

d{q,p) dx d{q,p) dy 



dH 

dx 

dH 
dy 

dH 

dz 



d{z,x) dt d{y,z) dt 

where we use the chain rule for the derivative of a function u{x, y, z) 

d _ , , du dx dudy du dz 

dt ^ ^ dx dt dy dt dz dt 



(A4) 



(A5) 



In case that d{x, y)/d{q, p) j^ 0, q and p can be expressed by x and y, and then the following 
equations are derived, 



dx dH d{x, y) 
dt dy d{q,p) 



dy 
dt 



dH d{x, y) 
dx d{q,p) 



(A6) 



using the relations 

d{q,p) 



d{x,y) \d{q,p) 



d{x,y) 



d{q,p) 

d{z,x) 







d{q,p) 
d{y,z) 



. 



(A7) 



In the same way, for the case of d{y, z)/d{q,p) y^ 0, q and p are expressed by y and z, and 
then the following equations are derived, 



dy dH d{y, z) 
'dt 



dz 



dH d{y,z) 



dz d{q,p) ' dt dy d{q,p) 

and for the case of d{z,x)/d{q,p) ^ 0, we obtain 



(A8) 



dx 



dH d{z,x) 



dz dHd{z,x) 

dt dx d{q,p) ' dt dz d{q,p) 



(A9) 
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Combing Eqs. (A6), (A8), and (A9), we can write down a set of equations, 

dx dH d{x,y) dHd{z,x) 

dt dy d{q,p) dz d{q,p) 

dy dHd{y,z) dH d{x,y) 

dt dz d{q,p) dx d{q,p) ' 

dz dH d{z, x) dH d(y, z) 

dt dx d{q,p) dy d{q,p) 

which is consistent with every expression of q and p. For example, in case that q = q{x, y) 
and p = p{x,y) {d{x,y)/d{q,p) 7^ 0), Eqs. (AlO)-(All) are reduced to Eq. (A6), and Eq. 
(A12) is equivalent to the Hamilton's equation of motion. 

Introducing a function G = G{x,y,z) that satisfies the conditions (12), Eqs. (A10)-(A12) 
are rewritten as the Nambu equations in the form of Eq. (4), 

dx d{H,G) dy d{H,G) dz d{H,G) 



(AlO) 
(All) 
(A12) 



dt d{y,z) ' dt d{z,x) ' dt d{x,y) 
These equations hold independently of the expression of q and p. 



(A13) 



B. Hidden Nambu systems in Nambu systems 

Let us formulate a Nambu mechanical system with an A^-plet Xj (i = 1, ■ ■ ■ , N) using N + r 
variables yj = yj{xi, ■ ■ ■ , xn) (j = 1, • " " > A^ + ^)) where r is a positive integer. We assume 
that at least r + 1 of {yj^, • • • , yj„}NB {jir " ^JN = ^r " ,^ + i~) do not vanish, where 
{yj^ , • • • , yjjv }nb is the Nambu bracket defined by Eq. (8). Then the equation for any function 
f{yi, ■ ■ ■ ,yN+r) = fixi, ■ ■ ■ ,xn) can be written as 

dj ^ d{f,H,,--- ,Hn-i) 
dt d{xi,X2,- ■ ■ ,xn) 

_ 1 v^ d{f,Hi,- ■ ■ ,Hn~i) f , , . 

where Haiyi,-" ,yN+r) = Haixi,--- ,xn) (a = 1, • • • ,N -I). 

Introducing functions Gc{yi, • • • , yN+r) = Gc{xi-, • • • , xn) (c = 1, • • • , r) that satisfy the 
relation 

1 v^ (9(Gi, ■ ■ ■ , Gr) _ J- -, , s 

„l 2^ ^jlJ2---JNJN + r--JN + r r,, \ — XV jl ■• V jll ' ' ' 1 V j N i^'B 1 \tiZ) 

' • • ,; _1 f^Vyiw + l) ■ • • li/jN + T-J 

jAf + 1,--- ,JiV + r — i 

Eq. (Bl) is rewritten in the form of Nambu equation, 

d~f d{f,Hi,--- ,Hn-i,Gu--- ,Gr) 



dt d{yi,y2r-- ,yN,yN+i,--- ,yN+ 
where we use the formula for any functions Aj = Aj {yi , • • • , yN+ 

d{Ai,A2,--- ,An,An+1,--- ,AN+r) 1 



(B3) 



N+r 



E 



d{yi,y2,--- ,yN,yN+ir-- ,yN+r) ^^-rl,. , , , 1 

^ d{Ai,A2, ■ ■ ■ , An) d{AN+i, ■■■ , An+t) ^^^. 
9{yji , yj2 , • • • , 2/j„ ) d{yj^^, , • • • , yj„^ J 
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By use of Eq. (B2), it can be shown that the Nambu bracket between Gc{xi, ■ ■ ■ ,X]\f) = 
Gcivi, ■■■ , VN+r) and arbitrary functions Ua{xi, ■■■ , xn) = Ua{yi, ■■■ , VN+r) (a = 1, • • • ,N 
1) vanishes such that 



N+r 



f^ ^ '^ Sr^ d{Gc,ui,--- ,UN~i) r , 



1 



N\ r\ 



N+r 

E 

d{Gc,ui,--- ,UN-i) d{Gi,--- ,Gr 



'Jlj2---JiVjN + l---JiV4 



X 



d{yj^ ,yj2,--- , vjm) d{yj n+1 ,■■■ , yj«+ J 

d{Gc,ui,- ■ ■ ,UN-l,Gi,- ■ ■ ,Gr] 



o[yi,y2, ■■■ ,yN, yN+i, ■■■ , yN+r) 

Hence Gc are constants. We can ehminate the constants by redefining Gc, and the resulting 
Gciyi,--- ,yN+r) = can be regarded as induced constraints, which are associated with 
enlarging the phase space from (xi, • • • , xn) to (yi, • • • , yj\f+r)- 

In this way, Nambu systems with N-plet xi (i = 1, • • • , N) can he formulated as Nambu 
systems with N + r-plet yj = yj{xi, • • • , X]\[) (j = 1, • • • , N + r). 
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